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ABSTRACT 

Evaporation  residues  produced  in  the  reactions  35Cl+169Tm 
and  37Cl+169Tm,175Lu,l8lTa  and  186W  were  observed  at  zero  degree 
utilizing  the  MIT-8NL  Recoil  Mass  Separator  (RMS).  The  recoiling 
nuclei  were  separated  from  the  beam  and  refocused  onto  a  surface 
barrier  detector  by  a  combination  of  electrostatic  and  magnetic 
fields  and  magnetic  quadrupole  lenses.  The  residual  nuclei  are 
alpha  radioactive  and  can  thus  be  identified  by  a  characteristic 
alpha  line  observed  after  the  arrival  pulse  of  the  evaporation 
residue.  The  recoiling  nuclei  also  passed  through  a  gas  ionization 
chamber  whereas  the  decay  alpha  did  not.  A  separate  anticoincidence 
spectrum  therefore  displayed  the  alphas  background  free.  Trends 
of  evaporation  residue  cross  section  were  charted  versus  Z2/A  and 
compared  to  statistical  evaporation  codes. 
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CHAPTER  I 

INTRODUCTION 

Since  the  Second  World  War,  the  alchemists  dream  of 
transmutation  has  become  a  commonplace,  if  not  very  economical  event. 
The  table  nuclides  has  been  increased  in  both  length  and  breadth  and 
there  are  now  more  than  a  thousand  known  nuclides.  The  table  has 
been  extended  far  past  the  fringe  of  the  line  of  stability  on  a 
curve  of  Z  vs.  N  to  include  very  large  nuclei  that  are  unstable  to 
spontaneous  fission  and  alpha  decay.  The  chart  has  also  been 
considerably  extended  toward  the  proton  drip  line  where  beta  decay 
predominates  on  the  fringe  and  alpha  decay  toward  the  center.  While 
there  is  no  expectation  of  finding  stable  nuclei  at  these  fringe 
areas,  their  production  serves  to  further  the  limited  body  of  knowledge 
concerning  the  systematics  of  the  stability  of  nuclei  to  various 
decays  as  well  as  further  refining  the  machinery  for  their  production 
and  detection.  The  other  direction  of  search,  at  higher  A,  has 
resulted  in  a  logjam  near  element  105  due  primarily  to  the  increased 
fission  competition.  The  various  methods  of  heavy  ion  fusion  to 
produce  these  more  and  more  massive  systems  have  resulted  in  little 
success.  This  may  be  the  result  of  theoretical  inadequacies  in  the 
fusion  reactions  producing  transfermium  nuclides  or  in  the 
methodology  and  apparatus  used  to  form  and  detect  them. 
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For  many  years  there  has  been  an  effort  to  produce  nuclides 
far  beyond  the  fringe  of  the  line  of  stability;  the  superheavy  nuclei. 
The  destabilizing  effects  of  the  large  Coulomb  repulsion  should  be 
overcome  by  the  stabilizing  effects  of  the  closed  proton  shell  at 
114,  116  or  possibly  126.  Many  avenues  of  approach  have  been  attempted 
to  achieve  fusion  coupled  with  survival  from  fission  in  this  area,  but 
to  no  avail  utilizing  any  of  the  present  heavy  ion  accelerators. 

One  method  is  to  bombard  a  very  heavy  nuclide  with  another 
very  heavy  one  such  as  uranium  on  uranium.  Here  it  is  hoped  that  the 
massive  compound  nucleus  would  fission  with  one  fragment  in  the  super 
heavy  region.  A  transfer  reaction  may  also  take  place  between  two 
massive  nuclei  producing  a  superheavy.  Some  modest  attempts  at  these 
reactions  have  been  unsuccessful  so  far  and  it  is  felt  that  the 
superheavy  fragment  may  be  produced  in  such  an  elongated  state  as  to 
fission  itself. 

Several  interesting  programs  have  been  undertaken  to  search 
for  super^eavy  nuclei.  The  Russians1  (Flerov,  1972)  have  studied 
medieval  stained  glass  windows  on  the  premise  that  a  superheavy 
nuclei,  chemically  similar  to  lead,  could  be  found  in  lead  glass.  A 
British  group  (Marinov,  et  al . ,  1971)  studied  tungsten  beam  stops 
from  the  30  GeV  CERN  proton  synchrotron.  They  estimated  that  a 
recoiling  tungsten  nucleus,  after  stopping  a  30  GeV  proton,  would 
have  sufficient  energy  to  fuse  with  another  tungsten  nucleus.  It 
was  hoped  that  a  long-lived  superheavy  nucleus  might  be  a  fission 
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product  of  such  a  reaction.  Terrestrial  ores,  lead  and  platinum,  have 
been  studied  and  nuclear  emulsions  exposed  to  high  altitude  cosmic 
rays  have  been  searched  for  extra-terrestrial  superheavies.  All  of 
these  attempts  have  failed. 

Another,  more  promising  method  has  been  to  bombard  a  heavy 
nuclide  with  a  medium  ion  such  as  *®Ca  (e.g.,  ^®Ca+^Cnr*296116*). 
Reactions  of  this  type  in  the  superheavy  compound  nucleus  region  have 
also  been  unsuccessful  so  far.  The  problem  has  been  the  dependence 
of  stability  against  fission  on  the  excitation  of  the  compound  nucleus. 
It  is  not  known  whether  production  of  a  •cold"  nucleus  to  minimize 
fission  is  favored  by  more  spherical  or  deformed  nuclei.  Currently 
a  modest  study  (below  uranium)  of  the  dependence  of  fusion  on  the 
excitation  of  the  compound  nucleus  is  being  conducted  by  the  M.I.T. 
Heavy  Ion  Group,  which  has  studied  the  reactions: 

65Cu  +  139la  -  204Rn  Egx  *  35-50  MeV 

35C1  +  169Tm  -  204Rn  Eex  =  42-72  MeV 

The  lack  of  success  in  producing  any  of  several  superheavy 

elements  is  cause  to  reevaluate  the  theoretical  and  experimental 

techniques  involved.  The  concept  of  the  stability  of  the  Coulomb 

barrier  may  need  revision  since  there  is  some  suggestion  that  under 

the  distortion  of  heavy  ion  interactions  the  increased  oblateness 

3 

of  the  projectile  results  in  an  increased  barrier.  The  liquid 
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drop  model  may  also  be  lacking  in  its  failure  to  predict  the  resistance 

4 

to  deformation  for  both  spherical  and  permanently  deformed  nuclei. 
Experimental  modifications  need  also  to  be  made  to  enhance  the 
detection  of  fusion  products  and  to  further  resolve  the  peculiar 
energies  of  their  decay. 

Compound  nucleus  formation,  using  a  relatively  light  projectile 
and  a  heavy  target  nucleus,  continues  to  give  copious  results  and 
that  is  the  method  followed  in  this  work  and  by  the  M.I.T.  Heavy  Ion 
Group . 

In  order  to  methodically  attack  the  problem  of  production  of 
superheavy  nuclei,  a  deliberate  approach  has  been  conducted  in 
compound  nucleus  formation  in  known  regions  of  the  chart  in  order  to 
gain  an  insight  into  the  mechanisms  of  compound  nucleus  formation  and 
to  validate  or  modify  current  theories. 

The  M.I.T.  Heavy  Ion  Group  has  aggressively  pursued  this 

(*] 

methodical  approach  and  has  conducted  or  has  approval  to  conduct  ' 
the  following  reactions: 
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These  reactions  have  been  selected  since  the  predicted  stable 
superheavy  elements  should  be  proton  rich.  With  the  exception  of 
the  Tungsten  Reaction  and  the  Dysprosium  group v  '  the  most  proton 
deficient  isotope  has  been  selected  to  give  a  proton  rich  compound 
nucleus  that  falls  in  the  area  of  alpha  decay. 

37  169  175 

The  experiments  described  in  this  work.  Cl  on  Tm,  Lu, 
181  186 

Ta  and  W  are  a  continuation  of  the  Heavy  Ion  Group's  program. 

This  set  of  reactions  was  chosen  for  its  gradual,  short  stepped,  march 

up  the  chart  of  nuclides.  This  work  is  an  attempt  to  quantify  the 
2 

Z  /A  dependence  of  the  fusion  cross  section  and  to  compare  it  with 
existing  theory,  primarily  the  ALICE  code.6  It  was  seen  immediately 
that  the  cross  section  diminished  rapidly  at  each  step  and  these  four 
reactions  were  taken  to  test  the  validity  of  the  ALICE  predictions. 

Here  is  where  the  limitations  of  the  detection  device  became 
paramount.  While  the  ALICE  code  predicts  a  certain  cross  section 
that  takes  into  account  all  compound  nucleus  evaporation  residues, 
the  "machine"  does  not  detect  them  all  due  to  its  limited  angular, 
charge  and  velocity  acceptance.  The  kinematics  program  BETH7  was 
utilized  to  specify  the  angula*'  and  velocity  distributions  of  compound 

O 

nucleus  formation.  The  code  RAYTRACE  was  used  to  quantify  the  charge 
and  velocity  acceptance  for  the  apparatus. 

The  term  "evaporation  residue",  above,  is  used  to  explain 
what  happens  to  the  compound  nucleus  after  fusion  and  coalescence  has 
occurred.  The  compound  nucleus,  highly  excited,  will  normally 


fission.  However,  a  small  fraction  will  survive  to  undergo  particle 
evaporation.  This  highly  excited  nucleus  then  de-excites  by  neutron, 
proton  and  alpha  particle  emission.  The  large  nucleus  left  after 
emission  of  a  few  particles  is  the  evaporation  residue.  These  then 
embed  themselves  in  a  solid  state  detector  and  decay  by  alpha  emission 
with  discrete,  identifiable,  energies.  All  experiments  are  conducted 
at  Brookhaven  National  Laboratory  utilizing  the  Tandem  Van  de  Graaff 
Generator  to  provide  a  copious  beam  of  ions  and  the  M.I.T.  Recoil  Mass 
Separator  (RMS)  to  separate  fusion  products  from  the  beam  particles. 
This  separation  is  necessary  since  simple  momentum  considerations  show 
the  fusion  products  to  be  forward  scattered  within  less  than  4°  on  the 
beam  line.  With  such  small  cross  sections  involved,  it  is  imperative 
to  move  these  precious  few  products  out  of  line  and  focus  them  into 
a  single  high  resolution,  solid  state  detector.  This  the  RMS  does. 

What  follows  is  but  a  small  portion  of  the  variety  of 
experiments  of  the  Heavy  Ion  Group.  Concurrent  and  complementary 
experiments  have  been  conducted  at  Berkeley  and  Darmstadt  and  it  is 
hoped  that,  as  the  body  of  knowledge  concerning  the  formation  and 
stability  of  compound  nuclei  increases,  and  as  larger  accelerators  are 
brought  on  line,  that  the  theory  of  superheavy  formation  will  come  to 
fruition. 


-18 


CHAPTER  II 

THEORY 


A.  General 

Any  group  of  theories  that  concerns  itself  with  such  a 
tremendously  complex  subject  as  heavy  nuclei  must,  of  course,  have  an 
empirical  basis.  Coefficients  and  form  factors  and  "normalizations" 
are  needed  for  a  theory  to  match  any  collected  data.  Some  theories, 
such  as  the  shell  model,  work  only  for  specific  nuclei  in  specific 
states.  Usually,  in  attempting  to  fit  theory  to  data,  or  vice  versa, 
more  questions  are  raised  than  answered.  This,  reasonably  enough, 
is  how  the  model  is  improved  and  is  a  basis  for  careful  experimentation, 
in  order  to  minimize  errors  in  the  refinement  of  existing  theory.  The 
work  presented  here,  hopefully,  both  adds  refinement  to  existing 
theoretical  models  as  well  as  improve  the  equipment  utilized  in  that 
refinement. 


B.  Reactions  Between  Complex  Nuclei 
When  two  complex  nuclei  enter  into  collisions,  two  main 
classes  of  phenomena  may  occur. 

When  the  system  follows  a  grazing  trajectory,  elastic 
scattering  is  the  dominant  effect  and  only  the  extreme  tails  of  the 
nuclear  matter  densities  overlap.  The  total  cross  section  is  a 
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function  of  the  summation  of  partial  waves  of  l  values.  Grazing 
interactions  are  related  to  the  largest  l  waves  whose  transmission 
coefficient  is  near  zero.  (See  Figure  1.) 

Figure  2  represents  the  two-body  interaction  potential 
restricted  to  one  dimension.  It  consists  of  the  nuclear  and  coulomb 
potentials  for  s-wave  interactions.  The  fusion  barrier  is  defined 
as  the  height  of  the  interaction  potential  at  zero  slope.  The 
intention  is  that  at  a  zero  slope  potential  the  system  can  be  trapped 
sufficiently  long  for  fusion  to  occur.  The  "effective"  radius  is 
defined  as  that  which  is  found  when  the  fusion  barrier  is  equated 
to  the  coulomb  potential.  Neglecting  the  nuclear  potential  will 
result  in  an  effective  radius  larger  than  the  actual  fusion  radius. 

As  target  and  projectile  charge  (and  therefore  coulomb  forces) 
increase,  a  larger  nuclear  potential  is  necessary  to  offset  the 
greater  Coulomb  repulsion. 

It  is  obvious  that  when  a  compound  nucleus  is  formed  a  large 

portion  of  the  kinetic  energy  is  transformed  into  excitation  energy 

of  the  compound  system.  This  fraction  is  equal  to  Aj/(A^  +  A2). 

Fragments  of  this  compound  system  possess  a  kinetic  energy  much 

smaller  than  the  entrance  channel  kinetic  energy,  showing  that 

g 

nuclear  matter  has  been  exchanged.  As  was  suggested  by  Cohen,  et 
al.  (1974),^  a  composite  nucleus  is  formed  for  a  system  that  "may 
well  correspond  to  an  intimate  contact  of  transient  amalgamation  of 
target  and  projectile."  If  this  fused  system  lasts  for  a  sufficiently 
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FIGURE'  2 


long  time,  it  is  called  a  compound  nucleus. 


Compound  nucleus  formation  has  several  features  that  allow  simple 

methods  of  analysis  and  yet  are  sensitive  to  particular  aspects  of  nuclear 

structure.  For  the  reaction  to  take  place,  the  projectile  nucleus  must 

2 

over  the  coulomb  barrier  between  itself  and  the  target.  E^g  =  Z-^e  /R 
where  R  s  is  the  separation  of  the  centers  of  the  two  ions. 

At  the  coulomb  barrier  energies,  the  corresponding  heavy  ion  wavelengths 
are  much  less  than  their  radii  so  that  their  motion  is  similar  to  that 
of  a  classical  particle.  (See  Table  I.) 


Absorption  of  a  complex  projectile  by  a  complex  target  means 
that,  after  a  time,  the  nucleus  of  both  take  on  new  characteristics 
in  a  single  nuclear  potential.11  The  total  mass  (A1  +  A 2)  recoils 
with  an  energy  ER  =  E.A^/(A^  +  A^)  in  the  direction  defined  by  the 
incoming  projectile  A-j  with  incident  energy  Ei .  Since  the  nucleus 
is  highly  excited  by  this  stage,  it  quickly  de-excites.  This  poses 
problems  for  its  detection  unless  the  path  of  its  de-excitation  is 
known  and  the  de-excited  "evaporation  residues"  are  easy  to  identify. 
From  the  simple  kinematics  above,  it  is  obvious  that  any  fused  nuclei 
will  be  forward  scattered  in  the  direction  of  the  beam  and,  of 
course,  have  a  much  greater  mass  (A^  +  A^)  than  either  the  beam  or 
any  knocked  out  target  nuclei.  The  angular  distribution  of  the 
compound  nucleus  is  within  quite  a  small  range  as  can  be  seen  in  the 
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TABLE  I 

Wavelengths  of  the  Projectiles  in  the  Laboratory  System 


Projectile 

(m) 

Target 

(M) 

Coulomb 
Barrier  E„ 
(MeV)  8 

Elab 

(MeV) 

*proj 

(fm) 

R  +  R  ** 
R1  R2 

(fm) 

*12c 

12c 

7.9 

15.7 

1.15 

6.59 

*40Ca 

,20Sn 

119.7 

159.6 

0.36 

12.03 

t  37 
d/Cl 

,69Tm- 

132.2 

161.1 

0.37 

12.76 

+  37C1 

186w 

138.9 

166.6 

0.36 

13.02 

*40Ca 

208Pb 

175.5 

209.2 

0.32 

13.46 

*208pb 

208Pb 

567.4 

1134.9 

0.14 

17.06 

♦Hodgson , 

Nuclear  Heavy 

Ion  Reactions, 

Ch.  4, 

Oxford  1978. 

^Present  Work 
**R.,,R2,  =  1.44  A1/3 
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angular  distribution  figures  accompanying  the  experimental  results. 
Nearly  all  of  the  distribution  is  within  8  degrees  of  the  beam  direction. 

Since  the  compound  nucleus  formation  cross  section,  o^,  is  less 
than  the  reaction  cross  section,  oR,  the  compound  nucleus  formation 
cross  section  is  therefore  the  sum  of  all  cross  section  up  to  l ^R: 

^CR 

°CN  ■  t\0  +  »>  h 


where  £CR  is  some  critical  angular  momentum  beyond  which  there  is  no 

more  contribution  to  compound  nucleus  formation  and  T „  is  the 

transmission  coefficient  for  a  particular  i  value.  One  of  the  simplest 

1 2 

modes  that  is  often  compared  to  experiment  is  due  to  Bass  (1973). 

The  critical  value  £CR  1S  determined  by  estimating  the  balance  between 

centrifugal ,  coulomb  and  nuclear  forces  in  the  two  body  system  at  the 

contact  distance.  Complete  fusion  is  presumed  to  occur  when  the 

attractive  nuclear  force  is  larger  than  the  sum  of  the  centrifugal  and 
1 3 

coulomb  repulsion.  At  this  point  the  total  potential  is: 


'TOT 


_  Z,Z2e‘ 


^CR^CR  +  1^2  ^SA/3AV3 


'12 


2yR 


12 


'12 


since  the  nuclear  potential  has  the  radial  form 

r  +  R 


a^'W3, 

VN(r)  *  - “o - -  exp 


'12 


12 


where  R12  *  i”q(a]^3  +  A^3) 
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as  is  the  surface  energy  constant  in  the  liquid  drop  mass  formula 

d  is  the  diffuseness  parameter 

y  is  the  reduced  mass 

(as  *  17  MeV,  rQ  =  1.07  fm,  d  =  1.35  fm) 


2 

This  relation  fits  the  increase  of  with  the  incident  energy 
until  saturation  is  reached 


* 


Z  2  R  a  _  7  7  dA^^A1^ 

■F  _  Z1Z2  .  49  R12asA1  A2  L\LZ*  asdA1  A2 
tSAT  ~  R12  +  25  2R^  R^ 

The  factor  49/25  comes  from  the  assumption  that  the  angular  momentum 
decreases  by  5/7  at  the  point  of  contact  (complete  sticking  and 
rigid  rotation  of  two  equal  masses  as  a  whole).  Above  no 

longer  depends  upon  energy  and  this  result  disagrees  with  experiment. 
The  sharp  conversion  of  radial  kinetic  energy  at  the  barrier  assumes 
that  there  is  no  friction  until  the  system  reaches  the  barrier  and 
then  infinite  friction  is  assumed.14  The  choice  of  rQ  3  1.07  is  also 
somewhat  arbitrary. 

Compound  nucleus  formation,  at  the  energies  and  masses  used 
in  this  study,  comprises  only  a  small  portion  of  the  total  reaction 
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cross  section.  Only  a  small  fraction  survives  fission  to  result  in 
evaporation  residues.  Due  to  the  deep  potential  well  to  charged 
particles  in  heavy  compound  nuclei,  the  primary  method  of  de-excitation 
is  neutron  emission  except  for  very  proton  rich  nuclei.  However  with 
increased  excitation,  the  compound  system  can  de-excite  more 
efficiently  through  the  emission  of  protons  and  especially  alphas  with 
their  high  binding  energy  and  large  angular  momenta.  Each  excitation 
function  rises  rapidly  from  zero  as  the  threshold  energy  is  passed, 
attains  a  broad  maximum  and  then  falls  due  to  competition  from 
reactions  emitting  more  neutrons.  (See  Figure  3.) 

The  energy  spectra  and  angular  distributions  of  particles 
evaporated  from  a  compound  nucleus  may  be  calculated  by  semiclassical 
statistical  theory  (Weisskopf,  1937)  providing  the  spin  of  the 
incident  particle  is  small  compared  with  that  of  the  compound 
nucleus.^  The  probability  of  emission  of  particle  y  with  energy 
and  spin  in  a  direction  9  is  given  by 


a(9,  Ey)  =  A 


2JTjdJ  X  (2SY+l)p(EXY,0) 


2LT^exp- 

o2+iT 

Jo 

fi  JL 

2a2 

2 

0 

L  yJ 

.  v, 

WJL(9Ey)dl 


where 


fiJl 

~T 

aY. 


-1 


K-D 

k 


(4k+l) 


'  ’ 

2 

r  ' 

(2k)! 

.  J  . 

iJL 

(2kk!)2 

2k 

ra 

P2kcos9 


and  Tj  and  T^  are  the  transmission  coefficients  in  the  entrance  and 
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c 


exit  channels,  J  is  the  spin  of  the  compound  nucleus,  L  is  the 

angular  momentum  of  the  emitted  particle  and  r j  is  the  total  decay 

width.  The  level  density  p(E  ,J)  of  the  residual  nucleus  at  the 

xy 

excitation  energy  E  and  spin  J  may  be  obtained  from  the  expression 
xy 

of  Long  (1966)^  and  Thomas  (1968)^  for  the  density  of  states  with 
energy  E  and  spin  J. 


p(E,J)  = 


12  at/4(En)='W)3/z  eXP 


2(ae) 


1/2 


*'9 


21 


2a 


where  t  is  the  nuclear  temperature,  a  is  the  single  particle  level 

2  2 

density  parameter,  a  -  Ift/h  the  spin  cut-off  parameter  and,  Ir 
the  rigid  body  moment  of  Inertia.  The  nuclear  temperature  is  related 
to  the  excitation  energy  E  by 


E  * 


+ 


mV 

21 


The  compound  nucleus,  in  its  excited  state,  emits  a  variety 
of  particles  (alphas,  neutrons,  protons  and  gammas).  However,  the 
probability  of  emission  differs  depending  on  both  the  angular  momentum 
and  excitation  energy  of  the  compound  system.  A  system  with  large 
angular  momentum  favors  alpha  particle  emission  since  each  alpha  can 
carry  off  more  angular  momentum  than  the  others.  When  the  excitation 
energy  of  the  compound  system  is  less  than  the  neutron  binding  energy 
(and  above  the  yrast  energy)  then  gamma  emission  is  favored. 

(Grover  and  Gelat,  1967)  Figure  4  shows  the  results  of  Grover's 
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calculations  for 


120 


Xe. 


C.  Statistical  Evaporation  Codes 

The  statistical  theory  of  Weisskopf  and  Ewing  is  utilized  by 
the  ALICE  code  to  describe  the  de-excitation  of  a  compound  nucleus 
through  particle  emission.  The  relation 


Pv(e)de 


(2s+l)v  uv  cav(e)p(Ef)de 


i  (2’+'>vk 

v*l 


v 


e  av(e)p(Ef)d£ 
'o 


describes  the  probability  of  emission  of  a  particle  v  with  energy 
between  e  and  e  +  de,  spin  S,  reduced  mass  y  and  inverse  cross 
section  a(e)  from  an  excited  nucleus  to  form  a  nucleus  of  level 
density  p{Ef).  This  density  can  be  approximated  by  the  Fermi  gas 

O  _ 

level  density  p(E,  0  ■  0)  *  E  exp  2/2aE  where  a  is  the  nuclear 

-1  * 

level  density  parameter  taken  to  be  about  A/10  MeV  .  The  assumption 
that  p(E,J)  =  (2J  +  l)p(E,  J  =  0)  is  also  made. 

The  evaporation  code  ALICE  calculates  the  probabilities  of 
emission  of  various  particles  for  each  possible  step  of  an  evaporation 
cascade.  This  then  leads  to  a  prediction  of  the  cross  section,  for 
a  specific  bombardment  energy,  for  the  eventual  evaporation  residue. 

It  also  calculates  the  total  reaction  cross  sections  from  the  inverse 


♦Calculations  for  this  study  varied  the  parameter  from  A/12  to  A/8. 


parabolic  model,  the  inverse  reaction  cross  sections  from  the  optical 
model  and  fission  competition  utilizing  the  rotating  liquid  drop  model 
of  Cohen,  Plasil  and  Swiatecki.^  As  noted  by  D.  Horn  (1976)^  the 
ALICE  code  fails  to  accurately  decrement  the  probability  distribution 
for  J  values  in  each  residual  nucleus  once  a  particle  has  been 
evaporated.  This  deficiency  is  partially  offset  by  a  program  option 
that  allows  specific  amounts  of  angular  momentum  to  be  removed  for 
each  kind  of  particle  emitted.  Utilizing  this  option,  D.  Horn 
determined  the  difference  to  be  only  minimal.  A  second  deficiency 
concerns  the  level  density  parameter.  The  level  density  was  calculated 
as  a  function  of  the  excitation  energy  E  rather  than  E  minus  the 
rotation  energy,  ErQt. 


I  *  Nuclear  Moment  of  Inertia 

This  more  exact  level  density  was  not  used  nor  the  deficiency 
calculated  in  this  study. 

The  inverse  cross  sections  were  calculated  using  the  optical 

model  with  a  volume  absorption  Woods-Saxon  potential,  which  utilized 

21 

the  global  parameters  of  Perey  (for  proton  evaporations)  and  of 
22 

Wllmore  and  Hodgson  (for  neutron  evaporation) .  Alpha  particle 

23 

parameters  were  from  the  potentials  of  Huizenga  and  Igo. 


* 
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The  following  parameters  were  used: 


A.  Neutrons 

V  »  47.01  -  0.267  E  -  0.00118E2  MeV 
WQ  *  9.52  -  0.053  MeV 
*  0  MeV 

rQ  «  1.322  -  (7.6  x  10~4)A  +  (4  x  10' V  -  (8  x  10'y)A'3  fm 
a  *  0.66  fm 

rQ  *  1.266  -  (3.7  x  10~4)A  +  (2  x  10'6)A2  -  (4  x  10'9)A3  fm 
a^  *  0.48  fm 

Wso  *  7-° 
rso  *  rD 
aso  *  aD 


B.  Protons 


a 

W0 

r0 

aQ 

Wso 

r$o 


=  53.3 

*  1.25 
=  0.65 

*  13.5 

*  1.25 
■  0.47 

*  7.5 
rQ,a 


0.55  E  +  27 (N-Z )/A  -  0.42  A1/3 


so 


=  ar 


C.  Alphas 


V 

w 

ro 

a 

r 


50.2 

12.3 

1.2  +  1.5  A 
0.564 

1.3 


-1/3 
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The  optical  model  assumes  that  the  nucleus  is  cranslucent  to  the 
incident  particle  wave.  The  projectile  therefore  has  a  specific  mean 
free  path  within  the  nucleus;  the  size  of  the  mean  free  path  then 
determines  if  absorption  occurs.  The  optical  potential  has  the  form: 


V  = 


Vc  -  Vf(x) 


■E 


Wf(x)w  - 


4WP 


£  F<y] 


h 

+ 

m  c 

w  J 

where  Vc  *  Coulomb  potential 
V  *  Real  potential 
W  =  Volume  absorption  potential 
WD  *  Surface  absorption  potential 
Vso  “  Spin  ort)it  potential 


f(x)  *  (1  *  eV 


*  (r 


rb> 


A1/3/a. 


As  previously  mentioned,  the  decomposition  of  the  total  cross  section 
into  partial  waves  is  given  by: 

«CN  *  »  JQ  <Z«  *  1  >TC 

where  the  transmission  coefficients  T^  are  related  to  the  barrier 
height,  B,  and  excitation  energy,  E,  by: 
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T£  »  [1  +  exp^CB-E)/^}]'1 


The  t»3 ^  are  the  vibrational  frequencies  of  a  harmonic  oscillator  having 
a  potential  energy  function  given  by  -V  and  are  related  to  the  real 
part  of  the  potential  by: 


1 

mdr2 


4 


where  V  is  the  total  potential  used  to  determine  the  total  fusion 
cross  section: 

where  VQ  =  67  MeV 

R  -  r0(A]/3  +  A^/3)  fm 
a  a  .57  fm 

The  ALICE  code  approximates  this  total  potential  by  an  inverse 

parabolic  potential  matched  in  position  and  curvature  at  its  maximum. 

Nuclear  binding  energies  are  calculated  by  the  use  of  the  Myers- 

24 

Swiatecki  Mass  formula. 

The  most  imperative  fit  that  must  be  made  with  the  ALICE 

code  Is  the  evaporation  to  fission  ratio,  or  the  total  evaporation 

25 

residue  cross  section.  There  are  over  45,000  combinations  of 
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options  that  can  be  selected  by  the  user  and,  given  a  set  of  potentials 

and  angular  momentum  parameters,  the  fission  barrier  must  be  adjusted. 

This  adjustment  is  made  by  varying  the  ratio  of  single  particle  level 

densities,  CLD,  the  level  density  parameter,  PLO,  and  the  barrier 

factor,  BARFAC,  in  the  ALICE  input  deck. 

The  most  important  parameters  of  the  statistical /rotating 

liquid  drop  model  are  the  ratios  of  the  single  particle  level  densities 

of  the  fissioning  nuclei  (at  the  saddle  point  deformation)  to  those 

of  the  residual  (evaporating)  nuclei  (at  the  equilibrium  deformation), 

26 

a^/av,  and  the  fission  barriers.  A  main  objective  of  this  work  is 

to  determine  a^/a^  ratios  and  fission  barriers  that  best  fit  the 

experimental  data.  Fission  barriers  B ( I )  are  given  by  the  differences 

between  rotating  energies  of  the  nucleus  at  the  saddle  point 

deformation,  Esp(I),  and  the  rotating  nucleus  at  the  equilibrium 

27 

deformation,  Emin(I). 

V"  *  W»  *B<'> 

The  barrier  factor  BARFAC  or  Bf  was  used  to  scale  the  rotating 
liquid  drop  fission  barriers,  B(I).  This  was  accomplished  by 
replacing  Esp(I)  by  Egp(I),  where 

EV'> '  W!)  * 


-36- 


Results  of  10  different  parameter  sets  are  shown  in  a 
later  chapter. 

Figure  5  shows  the  complete  fusion  cross  sections  for  the 
studied  reactions  as  generated  by  the  ALICE  code.  The  total  fusion 
cross  section  is  shown  only  to  establish  a  reference  for  1-value^, 
whereas  the  ratio  of  fission  to  evaooration  is  the  interesting 
quantity.  The  two  processes  are  in  competition  at  the  lower  1- values 
and  there  is  a  multiple  chance  for  the  compound  nucleus  to  fission 
after  each  successive  particle  boiloff.  Fission  is  completely 
dominant  at  the  higher  1-values.  Figure  6  shows  the  evaporation 
residue  cross  sections  for  the  studied  reactions.  Note  that  this 
represents  the  sum  of  all  of  the  possible  evaporation  channels  by 
which  the  compound  nucleus  can  de-excite.  Note  also  the  expanded 
scale  as  well  as  the  dashed  line  representing  the  average  complete 
fusion  cross  section  from  the  previous  figure.  For  these  reactions 
the  evaporation  residue  cross  sections  are  very  small  (total  cross 
sections  are  obta;ned  by  integrating  over  all  1-values)  as  compared 
to  the  fission  cross  sections. 

Table  II  shows  the  measured  and  calculated  evaporation 
residue  cross  sections,  the  latter  determined  by  the  ALICE  code 
for  10  different  parameter  sets. 

Having  determined  the  best  fit  set  of  entry  parameters  to 
obtain  a  fit,  or  bracket,  of  the  total  evaporation  residue  cross 
section,  a  check  must  be  made  to  determine  the  adequacy  of  the 


Angular  Momentum  (ih) 

FIGURE  6 


-40- 


individual  excitation  functions  for  each  evaporation  channel.  For 
this  study  the  lower  total  evaporation  residue  cross  section  limit  is 
given  by  parameter  set  =  0.80,  a^/av  *  1.03,  PLD  *  9.5,  will  be 
utilized.  This  set  was  chosen  as  the  closest  set  to  the  experimental 
data  that,  in  the  majority  of  the  cases,  was  less  than  the  experimental 
values.  Many  studies  have  shown  ALICE  to  overestimate  the  evaporation 
residue  cross  sections  and  it  is  hoped  that  this  study,  with  several 
low  values,  will  help  to  show  that  ALICE  can  be  an  effective,  and  cost 
effective,  tool  in  first  estimations  of  heavy  ion  fusion  reactions. 
Table  VIII  shows  the  calculated  and  measured  evaporation  residue  cross 
section  values  for  the  parameter  a^/a^  =  1.03.  The  table  includes  the 

OQ 

lighter  systems  (A^  »  98  to  172)  presented  by  Beckerman  and  Blann 
O  as  well  as  the  heavier  systems  from  this  work  (A^  =  204  -223).  The 

values  of  a^/av  =  1.03  continues  to  show  promise  in  predicting  total 
evaporation  residue  cross  section.  Note,  however,  that  the  fission 
barrier  factors,  B^.,  increase  with  A^.  This  constancy  of  the  ratio 
of  single  particle  level  densities  of  the  fissioning  nuclei  to  those 
of  the  residual  nuclei  as  well  as  the  linearity  of  the  fission  barrier 
factor  may  prove  helpful  in  predicting  the  survivability  of  heavier 
systems. 

Figure  6  shows  the  differential  cross  section  versus  angular 
momentum  shows  smooth  functions  for  all  reactions  except  the  Tungsten 
reaction.  This  compound  nucleus  differs  from  the  other  four  in  that 
it  posesses  an  odd  unpaired  proton.  The  Tungsten  target  is  an 
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even-even  nucleus  (resulting  in  an  odd-even  compound  nucleus)  while 

all  others  considered  were  odd-even  targets  (resulting  in  even-even 

compound  systems).  This  shell  effect  has  been  taken  into  account  by 

the  ALICE  code.  The  total  differential  cross  section  for  the  Tungsten 

reaction  is  not  smooth,  but  appears  with  two  intermediate  peaks.  These 

are  the  result  of  a  favored  evaporation  channel  (3na)  that  is  the 

major  contributor  to  the  total  evaporation  cross  section  and  that 

peaks  strongly  in  the  same  two  areas.  This  may  be  the  result  of  the 

223 

particular  deformation  of  the  compound  system.  Pa  (odd-even),  that 

favors  alpha  emission  at  specific  angular  momenta. 

The  assumption  made  by  ALICE,  that  the  emitting  nucleus  loses 

angular  momentum  at  the  rate  of  Ifi  for  neutrons  3ti  for  protons  and 

10h  for  alphas,  depends  very  much  on  the  size  of  the  compound  nucleus. 

At  high  A  this  assumption  becomes  even  less  valid.  The  code  MBIIZ 

handles  the  angular  momentum  loss  much  better  by  coupling  the  emitted 

* 

particle  angular  momenta  to  that  of  the  compound  system.  This  shows 
up  favorably  in  suppressing  the  multiple  charge  particle  boiloffs  that 
are  seldom  seen  (i.e.  3pXn  and  2aXn).  The  MBIIZ  code  is  costly  to 
run  and  should  be  used  only  as  a  refinement  upon  ALICE  data. 

0.  Velocity  and  Angular  Distributions  Generated  by  the 
SETH  Code 

Since  the  velocity  selector  of  the  RMS  passes  a  range  of 
velocities  rather  than  a  single  velocity,  an  experimental  distribution 


F(v),  is  a  convolution  of  the  true  velocity  profile,  g(v),  and  the 
velocity  dependent  solid  angle  function,  G(v-v'). 
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The  BETH  code  (Bolton,  1978)  calculates  the  velocity  distribution, 
g(v'),  of  the  particular  evaporation  residues.  Since  the  compound 
nucleus  de-excites  through  neutron,  proton  and  alpha  emission  (or  a 
combination  of  these  three),  the  evaporation  residue  (the  large 
nucleus  remaining)  recoils  in  a  distribution  about  the  center-of-mass 
velocity,  vcm.  Their  direction,  relative  to  the  beam  line,  is  also 
a  distribution  around  6=0°  rather  than  straight  along  it.  The  BETH 
program  assumes  that  particle  evaporation  is  sequential  and  isotropic. 
The  assumption  of  an  isotropic  distribution  of  evaporation  residues 
has  been  studied  by  Britt  and  Quintan"*0  who  found  that  the  anisotropy 
of  alpha  evaporation  is  small.  As  an  example  he  cites  that  the  ratio 
of  the  observed  anisotropy  for  10.5  MeV/nudeon  C  +  Au  in  the  center- 
of-mass  distribution  is  1.20  ±  0.20.  This  is  consistent  with 
evaporation  theory  presented  by  Ericson  to  first  order  in  ccs  0 


W(9)  =  1  +  (MaR2/I)(Erot/2T)cos20 


where  R  is  the  radius,  I  the  moment  of  inertia,  T  the  temperature  and 
Erot  the  rotational  energy  of  the  compound  nucleus.  Ericson 's 
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formula  gives  W(9)  as  approximately  1.3  to  1.4  (or  1.09  if  the  rigid 
body  moment  of  inertia  is  used).  The  BETH  assumption  of  isotropic 
emission,  then,  is  sufficient  to  give  a  general  overview  of  the 
angular  distribution  of  alpha  emission. 

The  following  parameters  were  used  by  C.  Bolton  in  the  BETH 

program: 


VELOCITY  DISTRIBUTION  IN  CM  FRAME 


d3a  ,  „  -k 

— r  (a,B,e,v) 

dvJ 


3  3 

^4  (a'.B'.e'.v')  (a",B" ,e",v")  d3v' 

dv,J  dvJ 


where  a  *  a'  +  a" 

(alpha  particles) 

S  «  8'  +  &" 

(protons ) 

e  =  £ 1  +  e" 

(neutrons) 

a'.S'.e’,  *  1  or  0 
a"  ,B" ,e" ,  *  1  or  0 
$  + 


INDIVIDUAL  (a,  p  or  nl  VELOCITY  DISTRIBUTIONS 


(001 )  *  Ee^E 


37  (010) 


r  _  -♦(E-E0) 

t/T  e 


E>Eo 

E<Eo 


rr 


-44- 


r 


37  0°°) 


-*<e-w  +  r 


ft  ftft7  e 


-YfE-Ej) 


E  < 

E  > 


Emax 

Emax 


where 


Y 

E 


1 


X 


Ek  -  5.2 

.961  Z/A1/3  (rQ  *  1.5  fm) 
.206  Z  +  2.53 
.9 


.24 
1/2(E 


max 


-  E.) 
-1 


Emax  +  (lnx) 

Emax  "  E1 
2<Emax  +  2  *  V 


4,r  were  determined  by  making  the  velocity  distribution  a  continuous 
fit  to  experimental  data. 

The  BETH  program  required  a  weighted  set  of  evaporation 
cascades  to  generate  the  function  q(v').  These  were  determined  from 
the  code  ALICE.  Once  generated,  each  distribution  is  added  (by 
weight)  to  yield  a  theoretical  velocity  profile  that  is  then 
transformed  to  the  lab  system  and  convoluted  with  n(v  -  v * )  to 
yield  a  velocity  distribution  that  can  be  compared  to  measured  values. 
The  transformation  is  the  Jacobian: 
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3  cos 8^  3V^ 

3  cose^  3V^ 

fVcJ 

3  cos0._  3V 
cm  cm 

3  V„  3  cose 
cm  cm 

w 

E.  Charge  Distribution 

To  determine  apparatus  settings  (covered  later),  an  average 

charge  state  q"  was  needed.  A  starting  point  for  determining  the 

distribution  was  to  use  the  formulae  of  Betz,  et  al .  and  Nikolaev 
32 

and  Dimitri ev 
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was  determined  as: 
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The  charge  state  determined  through  the  BETZ  formula  was  always 
found  to  be  too  low  by  three  to  fine  charge  states.  This  is 
consistent  with  the  BETZ  experimental  data  that  was  used  to  generate 
the  formula.  These  also  compare  with  the  RAYTRACE  acceptance  function 
shown  in  figure  12,  where  the  percent  change  in  charge  state  is  plotted 
versus  several  velocities. 


F.  Stability  of  Nuclei  as  a  Function  of  a  Critical  Value 
of  Z2/A 

Fission  may  be  viewed  as  a  study  of  nuclear  shape  deformations 

taken  to  an  extreme  limit,  resulting  in  separation  into  two  or  more 

fragments.  Only  binary  fission  will  be  discussed  here.  Adiabaticity , 

that  the  velocity  of  the  collective  motion,  here  the  fragment  of 

separation,  is  much  slower  than  the  individual  motion  of  the  nucleons, 

is  also  assumed.  The  use  of  the  liquid  drop  model  offers  a  first 

approximation  of  the  observed  features  of  fission  since  the 

33 

microscopic  motions  of  the  nucleons  are  not  accounted  for. 

The  surface  of  the  nuclear  drop  can  be  described  by: 


R(0)  =  R 


1  + 


N 

l 

n=2 


a  P„  cos  8 
n  n 


where  R  =  r.  A^  fm 

o 


The  coefficients  ap  are  generalized  deformation  coordinates 
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cons  trained  by  the  conditions  of  nuclear  incompressibility  and  the 
constancy  of  the  center-of-mass  coordinate  (origin) 


(Incompressibility) 
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The  surface  energy  is  proportional  to  the  surface  area  for 
a  spherical  drop: 


Es(6)  t  4ttR2  S 


where  S  is  the  surface  tension 


r r 


Es(a 


2  *a3  ’ 


. . )  =  4ttS 


R2(e) 


t  1  dRlell 
R(e)  de 


h  sinede* 


Jo 

if  the  deformation  parameters  are  small 


^a2’a3,#  ’ =  4*R0S 
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♦The  incompressibility  and  Constance  of  center  of  mass  equations 
require  ag  and  =  0. 
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The  coulomb  energy  of  a  charge  distribution  of  density  p(r) 
is 
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eGMtl  d;d;- 

|r- r'l 


with  a  uniform  charge  and  constant  volume  the  charge  density  (also 
constant)  becomes 


p  *  Ze/j  itR3  =  .414  Ze/A 


and 


Ec(°)  *  |  (Ze)2/1.2  A1/3 


For  small  distortions  the  coulomb  potential  then  becomes: 
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The  deformation  energy  is  the  combination  of  the  coulomb  and 
surface  terms. 
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By  restricting  the  potential  only  to  quadrupole  deformations  we 
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find 


AE  *  fcj  Es(0)0  -  x) 

where  x  is  the  fissility  parameter 

_  ZZ  3e2 
X'T40,n.2)JS 

For  x  <  1  the  drop  Is  stable  and  x  >  1  unstable  against  spontaneous 
fission. 


40  tt(1  .2) 
3e2 


40-50 


A  critical  value  of  Z  /A  becomes  the  dominant  factor  in  the 

stability  of  heavy  nuclei.  As  seen  in  figure  7  above  a  value  of 
2 

Z  /A  =  40,  the  spontaneous  fission  half  life  is  very  low.  There, 

the  probability  of  fission  versus  other  decay  modes  becomes  very 

high.  In  heavy  ion  fusion  reactions  in  the  regime  presented  in 
2 

this  work,  from  Z  /A  =  35.9  to  37.2,  the  nucleus  must  contend  with 
Its  own  natural  desire  to  fission  as  well  as  the  high  excitation 
energy  imparted  by  the  fusion  reaction  itself,  all  of  which  greatly 
increases  the  instability  of  the  compound  system.  This  results  in  a 
high  incidence  of  fission,  occurring  at  all  angular  momentum  values 
as  shown  by  ALICE  (see  figure  6},  with  very  little  remaining  to  be 
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de-excited  by  particle  emission.  This  is  the  point  of  departure  of 

2 

this  work;  to  quantify  the  stability  of  the  compound  systems  as  Z  /A 
increases  toward  a  critical  point.  While  that  critical  point  is  not 
reached,  the  trend  of  instability  is  quantified. 
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CHAPTER  III 

EXPERIMENTAL  APPARATUS 

A.  The  Recoil  Mass  Selector 

The  purpose  of  the  RMS  is  to  separate  fusion  products  from 
the  projectile  ion  beam.  The  beam  is  collimated  and  centered  onto  a 
thin  foil  target,  where  fusion  occurs.  The  majority  of  the  beam  is 
scattered.  This  scattering,  along  with  the  fusion  products  if  forward 
peaked  in  the  beam  direction,  however  the  fusion  products  have 
significantly  less  velocity  and  are  deflected  more  than  the  faster 
f-  moving  beam  by  a  series  of  electric  and  magnetic  fields.  These  fusion 

products  (plus  a  small  portion  of  the  beam  slowed  by  scattering)  are 
then  focused  onto  a  solid  state  detector  through  a  series  of 
quadrupole  magnets.  This  apparatus  is  well  suited  for  detection  of 
the  evaporation  residues  of  the  fusion  products  which  exhibit  a  small 
angular  distribution  and  a  low  yield.  Here,  zero  degree  separation 
is  a  necessity  and  is  accomplished  by  the  removal  of  the  evaporation 
residues  from  the  line  of  the  projectile  beam. 

The  beam  is  obtained  from  the  BNL  MP-7  Tandem  Van  de  Graaff 
Generator  with  a  maximum  terminal  voltage  of  about  14  MeV.  The  beam 
is  directed,  focused  and  shaped  magnetically  as  it  enters  the  target 
room.  A  collimator  consisting  of  two  aperture  sets  separated  by  a 
20  inch  pipe  is  the  first  portion  of  the  RMS  encountered.  The 
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apertures  are  mounted  on  a  base  that  allows  the  selection  of  a  variety 
of  sizes.  For  this  series  of  experiments  a  0.375  inch  diameter  was 
selected.  The  apertures  are  lined  with  gold  foil  to  minimize  any 
undesirable  reaction  products  that  may  result  from  the  beam  striking 
the  aperture  edge.  Electrical  connections  are  made  to  monitor  as 
current  any  beam  particles  impinging  on  the  collimator. 

The  collimated  beam  enters  the  target  chamber  (see  Figure  8) 

which  consists  of  a  round  drum  with  an  exit  port  cut  through  a  sliding 

o 

tape,  a  rectangular  slit,  a  target  ladder,  and  two  25  mm  solid  state 
detectors  (as  monitors  of  Rutherford  scattering  for  normalization 
purposes).  The  rectangular  slit  consists  of  four  independently 
adjustable  spectroscopy  grade  carbon  rods.  Current  is  monitored  from 
each  rod  and  optimal  beam  collimation  occurs  when  these  currents  are 
at  a  minimum.  Slit  dimensions  were  0.055  in.  vertical  and  0.150  in. 
horizontal  during  these  experiments.  The  target  chamber  has  several 
view  ports  to  inspect  the  interior  of  the  target  chamber  while  under 
vacuum.  The  target  ladder  is  mounted  one  inch  behind  this  final 
collimator  and  is  capable  of  supporting  four  target  frames.  This 
affords  the  ability  to  change  targets  without  interrupting  the 
vacuum. 

2 

Two  25  mm  solid  state  detectors  are  mounted  at  angles  of 
28.8°  and  distances  of  23.28  cm  from  the  beam  spot  on  the  target 
ladder.  This  results  in  a  solid  angle  of  acceptance  for  the  two 
monitor  detectors  of  0.07317  msr  (monitor  aperture  is  1/16  in.). 
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The  exit  port  of  the  target  chamber  (0.5  in.  diameter)  is 
mounted  on  a  sliding  steel  tape  pressed  against  a  neoprene  gasket  that 
is  sealed  to  preserve  the  vacuum.  This  equates  to  a  7.5  msr  solid 
angle  at  the  target  chamber  exit.  The  entire  recoil  mass  selector 
apparatus,  over  four  meters  long,  connects  to  this  two  inch  wide  tape 
and  1/2  inch  hole  and  can  rotate  about  the  target  spot.  This  enables 
the  RMS  to  be  utilized  for  a  variety  of  experiments  requiring 
observations  at  other  than  zero  degree. 

The  fusion  products  and  beam  exiting  the  target  chamber  first 
encounter  a  magnetic  auadrupole  pair  at  26  cm  from  target  (the  15/16 
in.  aperture*  here  equates  to  4.2  msr).  This  pair  of  quadrupoles  is 
designed  to  focus  the  exiting  particles  into  the  velocity  selector 
in  order  to  increase  the  acceptance  solid  angle. 

A  charged  particle  with  velocity  v  experiences  a  force  t 
when  moving  in  an  electromagnetic  according  to  the  Lorentz  equation: 

f  *  q (t  +  v  x  f) 


Thus  a  positively  charged  ion  moving  in  the  +z  direction 
will  experience  a  Lorentz  force: 


*n.b.  All  aperture  holes  mentioned  are  tapered  to  afford  easy  solid 
angle  calculations. 
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The  sign  of  the  force  component  is  significant  since  it  shows  a 
focusing  in  the  x  direction  and  a  defocusing  in  the  y  direction. 
Applying  the  thin  lens  approximation,*  valid  for  fg  »  l,  the 
effective  length,  the  focal  length  is  then 
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For  a  quadrupole  with  hyperbolic  pole  pieces. 
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where  R  is  the  aperture  radius  and  Bg  is  the  field  at  the  pole  tip. 

In  order  to  focus  a  beam  of  charged  particles,  two  quadrupoles 

are  utilized.  The  first  acts  a  converging  lens  in  the  x  direction 

* 

and  a  diverging  lens  in  the  y  direction.  The  second  acts  in  an 

A  A 

opposite  manner  diverging  x  and  converging  y  (see  Figure  9). 

The  pole  pieces  used  in  the  quadrupoles  are  6  in.  long  with 
an  aperture  diameter  of  3.0  in.  The  coils  are  water-cooled  and  each 
quad  is  protected  by  a  McDonnell  flow  switch  connected  to  the  power 
supply.  The  power  is  supplied  by  two  Hewlett  Package  6269B  power 
amplifiers.  The  purpose  of  the  first  pair  of  quadrupoles  is  to 
focus  a  larger  cone  of  particles  into  the  Wien  filter  than  would  be 

*n.b.  Thin  lens  approximation:  that  the  length  of  the  effective 

field  is  much  less  than  mv  /B„. 

o  q 
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avail  able  geometrically  over  the  distance  to  the  target.  The  Wien 
Fflter^’^  analyzes  the  charged  particles  according  to  velocity. 

The  fields  are  arranged  perpendicular  to  and  mutually  along  the  path 
of  the  particle.  The  Lorentz  force  will  vanish  if  its  velocity  is 
equal  to  the  ratio  of  the  Electric  to  the  magnetic  field  strengths. 
(vQ  *  E/B).  With  no  force  acting  on  it,  the  particle  will  not  be 
deflected.  On  the  other  hand,  a  particle  with  less  velocity  will  be 
deflected  up  while  the  one  of  higher  velocity  deflected  down.  This 
will  also  hold  true  in  the  fringing  fields  where  the  strengths  are  no 
longer  constant. 

The  fringing  fields  Introduce  complications,  however.  Upon 

encountering  the  velocity  selector,  the  particle  experiences  a 

magnetic  force  which  deflects  it  and  then  an  electric  force  to 

straighten  out  the  orbit.  However,  since  the  air  gap  is  larger  in 

the  magnet  than  in  the  electrostatic  deflector,  the  latter's 

* 

fringing  field  is  correspondingly  narrower.  While  the  electrostatic 
force  straightens  the  orbit,  it  offsets  the  trajectory  at  the 
entrance  and  exit  of  the  field.  This  offset  is  espressed  by 


where  Dm  and  are  the  magnetic  and  electric  air  gaps  and  I  is  the 
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fringing  field  shape  dependent  integral,  (see  Figure  10).  In 

the  present  device  I  =  0.85,  Dm  *  10  cm,  =  2.54  cm  and  B^/mvQ  *  169  cm 
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resulting  in  a  vertical  offset  of  0.47  cm  from  the  target  spot. 

There  are  two  sets  of  plates  in  the  velocity  selector.  First 
encountered  by  the  beam  is  an  electrostatic  field  which  deflects  the 
charged  particles  as  a  function  of  their  charge  to  energy  ratio 
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The  particles  are  allowed  to  drift  with  the  barely  deflected  beam 
particles  impinging  on  a  Faraday  cup  and  the  slower,  more  deflected 
evaporation  residues  entering  the  magnetic  field.  This  pure 
magnetic  field  has  the  same  effective  field  length  and  produces  the 
same  radius  of  curvature  as  the  electrostatic  field.  The  magnetic 
field  bends  particles  downward  (back  parallel  to  but  offset  from, 
the  beam  line)  according  to  the  charge  to  momentum  ratio, 
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The  third  portion  of  the  velocity  selector  is  the  t  x  § 
region  which  selects  velocities  such  that  vQ  =  E/B  is  undeflected. 

The  separation  of  pure  t  and  $  fields  by  the  drift  space 
results  in  a  net  upward  displacement  of  4.13  cm.  As  stated,  the 
more  energetic  beam  has  been  stopped  by  the  Faraday  cup.  Those 
particles  in  the  beam  tail  with  velocity  vQ  that  allow  it  to  pass 
through  the  filter  have  a  much  smaller  mass  and  are  deflected  into 
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an  adjustable  upper  baffle.  A  smaller  velocity  bite  will  result  with 
higher  fields.  Figure  11  shows  the  acceptance  solid  angle  versus  the 
percentage  of  velocity  deviation  from  vQ  for  p  =  169  cm  as  determined 
by  RAYTRACE. 

The  electrostatic  plates  are  supplied  by  Spellman  HVE 
Corporation  model  RHR  69P30/RVL  and  RHR  60N30/RVL  (±  60  KV).  The 
voltage  difference  is  remotely  monitored  using  a  digital  voltmeter. 

The  magnetic  fields  are  supplied  by  an  Alpha  (Donner  Systron)  440 
volt  supply.  A  digital  Hall  effect  gaussmeter  enables  remote  readings 
of  the  fields. 

Upon  leaving  the  velocity  selector  at  4.17  cm  above  the  beam 
line,  the  evaporation  residues  are  focused  onto  an  E-AE  telescope  by 
means  of  another  pair  of  quadrupoles.  The  pole  pieces  in  these 
quadrupoles  are  8  in.  long  and  the  aperture  diameter  is  4.125  in. 

These  are  also  water  cooled  and  flow  switch  protected.  The  ratio  of 
their  field  strengths  is  1.063  (downstream  to  upstream)  using 
resistors  in  a  current  dividing  circuit.  The  fields  are  monitored 
by  two  digital  voltmeters. 

The  Recoil  Mass  Selector,  to  recapitulate,  consists  of  a 
beam  collimator,  a  thin  target,  a  quadrupole  pair  to  focus  particles 
into  the  E  field  of  the  velocity  selector,  the  B  field  and  t  x  3 
region  of  the  velocity  selector  followed  by  another  quadrupole  pair 
that  focuses  the  separated  particles  onto  the  E-AE  detector.  The 
program  RAYTRACE  is  designed  to  show  the  path  of  a  specific  ray  of 


charged  particles  with  specific  charge  and  energy  that  is  emitted  at 
a  specific  angle  from  the  target.  RAYTRACE  takes  into  account  the 
abnormalities  of  the  fringing  field  areas  and  can  be  tailored  to  any 
configuration  of  multipole  magnets.  The  program,  then,  will  give  a 
velocity  bite  and  a  charge  bits  that  is  acceptable  to  a  specific 
configuration  of  multipoles  and  electromagnetic  fields.  Figure  12 
shows  a  contour  plot  Aq/q  vs.  Av/v  for  a  230  amu  particle  with  mean 
charge  state  q  =  17  and  energy  F  =  27  MeV.  Contour  lines  are  plotted 
from  10%  to  90%  acceptance  for  a  ±50  mr  ray  passing  through  the  system 
and  reaching  the  detector.  The  50%  acceptance  is  shown  to  be 
q  -  18%  to  q"  +  9%  and  v  -  5%  to  "v  +  3.5%.  Figure  13  shows  plots  of 
several  ±50  mr  rays  of  a  27  MeV  particle  to  show  the  action  of 
various  charge  states  and  velocities  in  the  RMS. 

To  determine  the  evaporation  residue  cross  section,  an 
angular  distribution  was  obtained  from  0°  to  6°  and  a  velocity 

i 

profile  measured  at  each  step.  The  velocity  distribution  is  measured 
by  varying  the  ratio  of  the  electric  to  magnetic  field.  All  of  the 
fields  are  varied  in  a  predetermined  fashioned  assuming  constant 
mass  but  a  mean  charge  state  that  varies  with  velocity.  Included 
into  the  velocity  distribution  is  the  velocity  acceptance  function 
of  the  instrument. 

For  momenta  and  charge  states  other  than  the  mean  charge 
and  momentum,  the  effective  solid  angle  of  acceptance  of  the  RMS 
rapidly  diminishes  because  of  chromatic  focusing  errors  (see 


FIGURE  12 
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Figure  IT).  The  effective  solid  angle  was  determined  by  elastically 
scattering  a  33  MeV  gold  beam  from  a  thin  gold  target  and  setting  the 
RMS  at  an  angle  equal  to  the  monitor  detectors  (28.8°).  The  ratio 
of  counting  rates  then  directly  gives  the  effective  solid  angle  of 
the  device  relative  to  the  solid  angle  of  the  monitors.  The 
elastically  scattered  gold  particles  have  a  sharp  velocity  spectrum 
but  the  charge  state  distribution  is  similar  to  that  of  the  evaporati 
residues.  Each  point  of  the  velocity  profile  represents  the  solid 
angle  integrated  over  the  charge  state  distribution. 

B.  The  Detector 

The  fusion  products  embed  themselves  in  a  solid  state 
detector  and  decay  there  primarily  by  alpha  emission.  If  the  alpha 
particles  are  stopped  in  the  depletion  layer  they  give  rise  to  full 
energy  peaks.  The  range  of  the  alpha  particle  is,  on  the  average, 

5.6  times  greater  than  that  of  the  heavier  evaporation  residue, 
and,  as  such,  43  percent  of  the  backward  emitted  alphas  escape 
through  the  surface  of  the  detector.  These  escaped  alphas  are 
registered  as  background  events  with  the  majority  of  the  energies 
far  below  the  lowest  energy  alpha  peak.  Therefore  a  correction  of 
1/.57  was  applied  to  cross  section  calculations. 

The  beam  rejection  capability  of  the  RMS  has  been  tested  at 
1 5 

1  part  in  10  particles  rejected  at  energies  above  30  MeV.  In  the 
6-30  MeV  range  a  few  background  events,  possibly  due  to  be&m 
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particles,  were  recorded.  Low  energy  scattered  beam  particles 
(approximately  200  sec"^ )  with  a  velocity  comparable  to  that  of  the 
fusion  products  (energy  equal  to  5  MeV),  however,  do  pass  through  the 
system. 

Since  the  low  energy  beam  tail  interferes  with  the  low  energy 
portion  of  the  alpha  spectrum,  a  AE  gas  ionization  chamber  is  placed 
in  front  of  the  E  detector.  The  AE  detector  is  used  to  create  an 
anticoincidence  signal  that  eliminates  the  background  since  only  an 
E  signal  resulting  from  a  decay  wi thin  the  detector  rather  than  an 
E  signal  due  to  an  entering  particle  is  desired.  Figure  14  shows 
two  spectra,  the  too  without  and  the  bottom  with,  the  use  of  the  AE 

chamber  in  anticoincidence  with  the  E  counter. 
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The  AE-E  detector  is  a  machined  aluminum  gas  chamber  with 

a  detector  holder.  The  surface  barrier  detector  is  held  in  place  by 

a  rear  mounted  Microdot  connector.  The  front  face  of  the  detector 

serves  as  the  rear  cathode  of  the  proportional  counter  whose  anode 

wire  is  mounted  just  outside  the  active  volume.  The  10  micron 

nichrome  wire  is  glued  to  the  contact  posts  under  5.2g  tension  and 

the  glue  points  painted  with  a  silver  conductor  to  insure  proper 

2 

contact.  The  window  consists  of  60-100  ug/cm  polypropylene  on  a 
removable  frame  that  fits  onto  a  special  cover  plate  for  the 
counter  (see  Figure  15).  The  Markham  counter,  here  used  only  to 
generate  a  AE  pulse,  can  also  be  utilized  as  a  two  dimensional 
position  sensitive  detector  as  outlined  in  the  referenced  article. 


e-Ae  counter 
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The  E  detector  was  calibrated  by  an  Americium-241  alpha 
source.  Energy  shifts  due  to  dead  layers  on  the  detector  and  source 
surface  were  studied  by  orienting  each  separately  at  60°  to  the 
normal  of  the  other.  This  produced  a  30  keV  shift  in  the  peak  energy. 
The  system  electronics  were  tested  for  linearity  with  a  pulser  to 
determine  slope  and  zero  point  energy.  Since  the  alpha  decay  occurs 
in  the  depletion  layer,  the  recoil  of  the  heavy  nucleus  that  decays 
also  contributes  to  the  pulse  height.  Approximately  12%  of  the  recoil 
energy  (-130  keV)  is  registered  however. 

C.  The  Electronics 

The  input  into  the  BNL  E7  computer  consists  of  the  AE-E 

system  and  a  timing  system.  Figure  16  shows  the  AE-E  system.  An 

"E"  signal  and  an  "All"  signal  are  displayed  as  well  as  the 

anticoincidenced  alpha  gate  in  order  to  be  able  to  see  an  overview 

» 

of  the  total  spectrum  as  well  as  to  see  where  the  alpha  spectrum 
falls  between  the  beam  tail  and  the  knockouts.  The  delay  amplifier 
on  the  E  signal  is  used  to  adjust  the  relative  timing  between  the 
E  pulse  and  the  gate  pulse.  The  Gate  and  Delay  Generator  (GDG) 
along  the  AE  signal  is  used  to  broaden  the  square  peak  put  out  by 
the  Timing  Single  Channel  Analyzer  (TSCA)  to  4.5  usee,  so  that  any 
system  fluctuations  are  not  allowed  to  disrupt  a  coincident  signal. 

The  final  GDG  reduces  the  pulse,  after  passing  through  the  universal 
coincidence  box,  to  0.5  ps. 


TIMING  ELECTRONICS 
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The  other  portion  of  the  electronics  is  the  timing  mechanism 
(Fig.  17).  Here,  the  evaporation  residue  pulse  was  used  to 
simultaneously  start  a  mul tf scalar  and  trigger  an  electrostatic 
deflector  that  cut  off  the  beam.  Data  is  recorded  in  a  2-D  time  vs. 
energy  mode.  The  beam  is  held  off  for  a  duration  equal  to  several 
half  lives  of  the  parent  nuclei.  The  2-D  timing  spectrum  can  be 
gated  to  up  to  16  separate  energy  peaks  by  the  E-7  computer.  This 
allows  for  the  simultaneous  timing  of  16  separate  alpha  spectra  which 
are  very  useful  in  determining  the  decaying  nuclide  when  similar 
alpha  energies  are  present.  An  accurate  half  life  is  also  necessary 
in  determining  isomerism. 
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CHAPTER  IV 


EXPERIMENTAL  RESULTS 


A.  37C1  ♦  169Tm 

This  experiment,  utilizing  the  very  plentiful  37C1  beam  (at 

the  BNL  Tandem),  is  included  in  this  work  as  a  point  of  departure.  The 

author  joined  the  research  group  after  much  of  the  experimental  work 

had  been  accomplished  and  participated  only  in  obtaining  velocity 

distributions  and  lifetime  data  and  performing  the  various  ALICE 

calculations.  It  is  included  in  this  work  since  all  system  improvements 

were  made  after  the  primary  data  was  obtained  and  because  all 

calculations  in  this  work  (to  obtain  absolute  cross  sections)  are  based 

on  method  conceived  for  the  Thulium  experiment.  The  soon  to  be 
30 

published  paper  detailing  this  experiment  contains  both  the 
results  and  the  methodology  of  the  experiment  and  will  only  be 
partially  reproduced  here. 

To  determine  an  absolute  cross  section  for  the  production  of 

a  compound  nucleus  that  survives  fission,  an  integration  must  be 

accomplished  over  both  an  angular  and  a  velocity  distribution.  It 

was  shown  in  figure  11  that  the  total  velocity  accpetance  of  the  RMS 

ran  from  +6%  to  -6%  of  the  expected  velocity.  This  was  obtained  by 

the  elastic  scattering  of  a  gold  beam  matched  in  energy  to  the 

2 

fusion  products  at  an  angle  equal  to  that  of  the  50  mm  monitor 


-75- 


detectors  used  for  normalization.  This  established  an  effective  solid 
angle  of  acceptance  of  the  RMS.  It  was  also  necessary  to  determine 
the  average  detector  to  monitor  counts.  This  was  accomplished  by 
superimposing  figure  11  over  the  several  intervals  of  the  velocity 
profiles  and  obtaining  an  average  value.  This  is  shown  in  figure  18 
for  a  3  point  velocity  profile.  The  integration  is  then  accomplished 
by  approximating  the  integral 


a  =  2tt 
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where  f  is  the  scaling  factor  obtained  from  the  normalized  count  to 
the  solid  angle  of  acceptance  ratio  and  A0  =  2°. 

In  the  conduct  of  an  experiment  to  study  the  various 
evaporation  residues  of  a  particular  compound  system,  a  separate 
velocity  spectrum  is  needed  for  each  particular  species  of  evaporation. 
It  is  impractical,  however,  in  a  study  of  several  reactions,  each  at 
several  Ub  energies,  to  obtain  a  detailed  velocity  and  angular 
distribution  at  every  energy.  In  the  conduct  of  these  experiments, 
one  velocity  distribution  and  one  angular  distribution  was  obtained 
for  each  evaporation  residue  at  a  lab  energy  near  the  peak  in  the 
total  cross  section.  For  the  ^C1  +  ^Tm  reaction  a  detailed 
velocity  profile  of  12  points  was  taken  ranging  from  -14. 7%  to 
+20.7%  of  vQ.  In  all  other  reactions  a  modified  3  point  velocity 
profile  of  +10%,  0,  -10%  of  vQ  was  taken.  All  angular  distributions 
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were  taken  in  2  degree  steps  from  0  to  8  degrees.  The  absolute  cross 
sections  for  each  evaporation  residue  were  then  obtained  by 
integrating  the  velocity  bite  taken  at  each  angle  over  the  8  degrees 
of  the  angular  distribution.  This  resulted  in  an  absolute  cross 
section  for  each  evaporation  residue  at  one  particular  lab  energy. 
Then,  based  on  the  ratio  of  normalized  counts  to  cross  section  at 
that  energy,  absolute  cross  section  values  could  be  scaled  to  other 
lab  energies.  Statistically  this  method  is  satisfactory  when  the 
angular  and  velocity  distributions  are  taken  at  the  peak  of  a 
particular  evaporation  residue  excitation  function;  however  when  the 
cross  section  must  be  scaled  from  the  tail  of  a  distribution,  a 
correspondingly  large  statistical  error  is  carried  forward.  This  is 
displayed  in  the  large  errors  in  some  of  the  low  cross  section 
reactions . 

Figure  19  shows  the  velocity  distributions  for  a  variety  of 
evaporation  residues.  The  compound  nucleus  with  a  forward  velocity 
in  the  CM  frame,  will  change  velocity  with  each  evaporation.  A 
pure  neutron  cascade  (i.e.  4n)  will  produce  little  deviation  in 
velocity  since  there  is  no  coulombic  barrier  to  overcome.  The 
charged  protons  and  alphas  impart  a  much  greater  change  in  velocity 
and  their  respective  distributions  are  therefore  not  center  peaked. 

A  difficulty  arises  in  determining  whether  a  particular 
evaporation  residue  is  a  daughter  (alpha  decay  in  the  detector)  or 
the  result  of  a  prompt  alpha  emission.  The  evidence  for  the  prompt 
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alpha  is  obtained  by  subtracting  the  parent  distribution  (i.e.  4n) 
from  the  daughter  (i.e.  4ncx);  the  resultant  distribution  (not  center 
peaked)  is  that  of  the  prompt  alpha  emission.  This  is  shown  in 
figure  19  where  the  6.636  MeV  (upper  left)  parent  is  subtracted  from 
the  6.180  MeV  daughter  (lower  left)  resulting  in  the  4nla  parent 
evaporation  residue  (dotted  line).  All  prompt  alpha  evaporation 
residues  in  this  study  were  obtained  in  this  manner. 

Figure  20  shows  the  angular  distributions  for  several 
evaporation  residues.  Here  also  the  prompt  alpha  angular 
distributions  were  determined  in  the  same  manner  as  the  velocity 
profiles.  Figure  21  shows  the  BETH  program  analysis  for  comparison. 
It  should  be  noted  that  the  BETH  results  are  an  integration  over  a 
wider  than  normal  velocity  profile.  This  results,  especi*  !y,  in  the 
alpha  case,  in  a  broadening  of  the  angular  distribution  and  an 
increase  in  the  zero  degree  valpes.  BETH  calculates  the  distribution 
only  on  a  kinematic  basis  with  no  regard  for  nuclear  or  RMS  effects. 
Measured  results  are  less  distinct  due  to  energy  loss,  straggling 
and  multiple  scattering  in  the  target. 

Figure  22  shows  the  excitation  functions  for  the  parent  and 

daughter  evaporation  residues  plotted  versus  lab  energy.  Note  the 

"tracking"  of  each  daughter  under  its  respective  parent.  The 

202 

predominant  peak  at  Rn  (4n)  is  the  major  contributor  to  the  total 
cross  section  plotted  in  figure  23  with  a  peak  evaporation  cross 
section  at  1600  pb  ±  300  at  165  MeV. 
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B.  35C1  +  169Tm 

2 

This  reaction  was  selected  as  a  midpoint  reaction  in  the  Z  /A 
distribution  and  to  show  the  effects  of  a  small  decrease  in  the  neutron 
number  (and  hence  an  increase  in  the  value  of  Z  /A).  This  reaction 
was  also  included  in  an  analysis  of  the  effects  of  excitation  energy 
and  angular  momentum  on  evaporation  residue  formation.  In  that  study 
the  same  compound  nucleus  (^Rn)  was  formed  through  the  reactions  of 
®^Cu  +  ^La  and  33C1  +  ^Tm.  The  copper  reaction  is  not  included 
in  this  study  since  the  data  is  incomplete.  A  much  more  complete  set 
of  excitation  functions  was  obtained  for  the  chlorine  reaction. 


Figures  25  to  27  show  the  angular  distributions  for  several 
of  the  parent  evaporation  residues.  Figure  25  shows  the  contrast 
between  the  pure  neutron  and  the  neutron  and  alpha  evaporation 
angular  distributions.  The  large  values  at  angles  4  and  6  degrees 
are  the  major  contributors  to  the  high  cross  section  values  for  the 
alpha  evaporation  residues  since  the  multiplication  of  the 
differential  cross  section  by  sin  9  near  zero  degree  results  in  a 
near  zero  contribution  to  the  absolute  cross  section.  (See  Table  III.) 
The  velocity  distribution  for  this  reaction  was  obtained  us*'ng  the 
three  point  method.  The  value  of  ±10%  of  vQ  was  chosen  since  the 
peak  of  alpha  evaporations  should  occur  in  that  region. 

Figures  28  to  31  show  the  excitation  functions  of  the 
several  evaporation  residues  of  this  reaction.  Figure  28  is 
significant  in  that  it  is  very  similar  to  the  3^C1  +  ^Tm  reaction 
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in  both  the  position  of  the  excitation  functions  and  their  relative 
peak  heights.  Of  further  interest  are  the  sets  of  other  excitation 
functions  involving  alphas  and  protons  mixed  with  neutrons  in 
evaporation.  Here,  as  in  several  other  reactions,  the  4np  and  the 
3na  reactions  are  the  strongest.  Also  of  note  is  the  presence  of  the 
2n2p  and  3n2p  evaporation  modes.  These  are  purposely  listed  as  shown 
rather  in  an  alpha  designation  since  they  occur  at  high  excitation 
energies  and  in  regions  where  4  or  5  particle  evaporations  occur. 

(See  Figures  29  to  31.)  The  total  evaporation  residue  cross  section 
in  shown  in  Figure  32.  The  large  errors  associated  with  the 
160-165  MeV  lab  energy  range  are  due  primarily  to  the  4n  excitation 
function  since,  as  previously  mentioned,  the  cross  section  scaling 
was  based  on  a  low  4n  value  at  a  lab  energy  of  175  MeV  and  the  error 
was  correspondingly  scaled.  In  Table  III  is  shown  the  individual 
cross  sections  that  contribute  to  the  total  evaporation  residue  cross 
section  function. 

199 

This  reaction  produced  the  new  isotope  Rn  and  its 
199n  40 

metastable  state  Rn.  It  was  possible  to  predict  the  expected 
alpha  decay  energy  by  charting  the  decay  ssystematics  of  the  neighboring 
Radan,  Astatine,  Polonium  and  Francium  isotopes.  Figure  24  shows 
the  decay  systematics  and  the  position  of  the  ^99Rn  and  ^"mRn 
isotopes  (open  circles).  Also  shown  are  the  alpha  decay  energies  of 
the  Polonium  daughters.  New  alpha  decay  energies  are: 


Rn  6.99  MeV 
199mRn  7.06  MeV 


”CI  +  '°*Tm  -  204  Rn 
Angular  Distribution 


•  6.380  3n  a 
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FIGURE  25 
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35CI  +  l69Tm  -  204 Rn 


Angular  Distribution 
•  6.636  4np 
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FIGURE  26 


Relative  Counts 
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FIGURE  27 
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TABLE  III 

Absolute  Cross  Sections  (in  ub)  for  ^C1  +  ^Tm 


Wtev) 

E0(HeV) 

150 

155 

160 

165 

170 

175 

180 

185 

2n2p 

9.4 

64.3 

98.1 

31.1 

3.1 

5.859 

±3.1 

±21.4 

±32.7 

| 

±10.4 

±1.0 

3n2p 

4.4 

10.4 

25.6 

44.0 

29.2 

6.063 

±1.2 

±2.9 

±7.1 

±12.2 

±8.2 

2np 

13.9 

58.8 

45.0 

28.9 

10.9 

2.9 

1.3  i 

6.354 

±9.8 

±41.5 

±31.8 

±20.4 

±7.7 

±2.1 

±.9 

3np 

5.7 

23.2 

100.0 

74.4 

47.1 

17.2 

4.2 

6.415 

±.  7 

±3.0 

±12.9 

±9.6 

±6.1 

±2.2 

±1.8 

3np 

6.2 

16.9 

11.2 

7.1 

4.2 

6.460 

±2.5 

±6.9 

±4.6 

±2.9 

±1.7 

4np 

15.6 

77.5 

155.5 

175.9 

77.8 

6.636 

±.9 

±4.4 

±8.7 

±9.9 

±4.4 

3n 

7.3 

53.8 

17.2 

6.77 

±4.2 

±31.1 

±9.9 

4n 

20.0 

93.1 

111. 9 

20.3 

3.2 

1.8 

6.900 

±11.5 

±53.8 

±64.6 

±11.7 

±1.8 

±1.0 

5n 

.5 

2.0 

.5 

6.990 

±.l 

±.  7 

±. ! 

5n 

.5 

2.4 

4.4 

3.9 

.5 

7.060 

±.2 

±.8 

±1.0 

±1.4 

±.2 

3na 

28.4 

105.9 

181.1 

145.1 

44.0 

6.380 

±5.5 

±20.4 

±34.9 

±27.9 

±8.5 

°T0T 


21.1 

±14 


138.3  222.5  448.5  486.8  424.0  295.9  111.9 
±84.8  f±109.8  {±148.9  (±109.4  ±68.8  ±38.9  ±14.6 
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C.  37C1  +  175Lu 

This  reaction  begins  the  series  of  reactions  with  very  low 

37 

cross  sections.  The  very  plentiful  Cl  beam  was  used  in  order  to 

obtain  beam  currents  of  200  na.  The  high  beam  current  was  chosen  to 

optimize  the  fusion  count  rate  even  though  it  resulted  in  a  nucleus 

41 

with  a  high  excitation.  One  study  of  the  program  ALICE  showed  that 

high  excitation  energy  resulted  in  a  decreased  evaporation  residue 

production.  The  recommended  projectile  nucleus  to  minimize  excitation 

A  ?  a 

CN  6 nCN 

energy  was  in  the  range  of  -j-  or  — j— »  however  the  intensity  factor 
of  100  gained  by  using  the  200  na  beam  greatly  offset  the  better 
production  efficiency  (larger  cross  section  a^R)  of  the  heavier  beam. 

Figures  33  to  35  show  the  angular  distributions  of  the 
various  evaporation  residue  modes.  In  the  regime  of  this  reaction, 
pairs  of  the  evaporation  residues  by  accident,  exhibit  identical  alpha 
decay  energies.  Separate  evaporation  modes  (i.e.  4n  from  5n)  could 
not  be  distinguished  by  using  velocity  or  angular  distributions  as  in 
the  prompt  alpha  case,  but  through  the  energy  (lab)  distribution  of 
their  excitation  functions.  This  is  shown  by  the  dashed  lines  in 
figures  36  to  39.  Angular  distributions  again  exhibit  the  characteristic 
shapes  of  their  respective  species.  This  constancy  of  distribution 
assists  in  the  Identification  of  individual  evaporation  residues  and 
In  the  certification  of  prompt  alpha  evaporation  modes.  Figure  40 
shows  the  total  evaporation  residue  cross  section  for  the  37C1  +  ^7\u 
reaction.  Individual  cross  sections  are  shown  in  Table  IV  for  each 
mode  of  evaporation. 


Relative  Counts 
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FIGURE  33 
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FIGURE  38 


FIGURE  40 
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TABLE  IV 

Absolute  Cross  Sections  (in  pb)  for  ^C1  +  ^tu 


Eub(nev) 

Ea(HeV) 

150 

155 

160 

165 

170 

175 

180 

185 

190 

3n2p/2n2p 

.5 

1.8 

8.5 

11.5 

21.5 

3.5 

9.8 

9.0 

6.129 

±.l 

±.3 

±2.3 

±3.0 

±5.7 

±.8 

±2.5 

±2.3 

2np/3np 

4.3 

8.5 

15.2 

33.1 

31.4 

12.4 

8.1 

2.2 

6.635 

±.5 

±1.1 

±1.9 

±4.2 

±4.1 

±1.6 

±1.1 

±.3 

4np 

2.2 

24.4 

35.3 

41.1 

29.7 

6.159 

±.2 

±2.5 

±3.5 

±4.1 

±3.0 

3n 

.72 

37.7 

62.2 

27.1 

11.5 

2.4 

.9 

7.002 

±.6 

±16.4 

±27.0 

±11.8 

±5.0 

±1.0 

±.4  j 

4n/5n 

1.0 

20.3 

48.0 

43.7 

27.8 

15.8  I 

7.8 

.3 

7.126 

±.l 

±1.6 

±3.9 

±3.6 

±2.3 

±1.3 

±.6 

±.2 

6n 

.6 

.8 

.8 

.2 

7.314 

±.2 

±.3. 

±.3 

±.l 

3na* 

34.9 

14.8 

6.261 

±10.3 

±4.4 

4na 

14.7 

34.2 

32.** 

16.4 

6.436 

±3.6 

±8.4 

±8.1 

±4.0 

-? 
o 
— i 

.72 

43.5 

92.5 

99.1 

136.9 

137.0 

102.2 

99.6 

60.3 

±.6 

±17.1 

±30.0 

±19.9 

±26.3 

±23.6 

±16.0 

±16.4 

±9.8 

♦The  reported  cross  section  is  10%  of  the  calculated  midpoint.  This 
was  done  due  to  the  large  error  introduced  in  integrating  over  the 
angular  distribution. 
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0.  37C1  +  181Ta  and  37C1  +  i86W 


These  two  reactions  are  considered  together  since  their  Z  /A 
values  are  very  similar  and  since  the  tungsten  data  is  less  than 
complete.  Once  again  in  this  region  several  of  the  evaporation  residues 
had  similar  alpha  decay  energies.  These  are  noted  by  dashed  lines  in 
figures  44  to  46.  Tungsten  excitation  functions  appear  in  figures  50 
and  51.  Angular  distributions  for  the  tantalum  appear  in  figures  41 
to  43.  Angular  distributions  for  the  tungsten  data  scaled  to  the 
37C1  +  ^Tm  shape,  appear  in  figures  48  and  49.  Tables  V  and  VI  show 
the  individual  cross  sections  for  the  two  reactions  and  total 


evaporation  residue  functions  appear  ir.  figures  47  and  52. 

Both  the  tantalum  and  the  tungsten  reactions  were  of  such  low 

cross  section  that  the  detection  system  was  severely  tested.  Only 

through  the  use  of  the  anticoincidence  counter  could  evaporation 

residues  with  such  low  cross  sections  be  observed  (in  some  cases  as 

♦ 

low  as  12  nanobarns).  Alpha  evaporations  begin  to  dominate  in  this 
region.  The  increased  excitation  of  the  nucleus  is  more  easily 
dissipated  through  the  evaporation  of  an  alpha  particle  rather  than  a 
neutron.  The  tungsten  reaction  displayed  many  problems  that  will 
surely  be  present  in  future  experiments  at  higher  Z.  In  many  instances 
the  half  lives  for  alpha  decay  of  the  evaporation  residues  were  of  the 
order  of  one  microsecond.  This  allowed  the  compound  nucleus  to 


evaporate  and  then  for  the  evaporation  residue  to  decay  before  striking 
the  detector.  Oue  to  the  time  constraints  and  a  low  count  rate,  no 
angular  distribution  was  obtained  in  the  tungsten  reaction.  To 


facilitate  the  determination  of  an  absolute  cross  section,  the  angular 
distributions  of  similar  evaporation  residues  in  the  ^C1  +  ^Tm 
reaction  were  used.  Excitation  functions  of  the  short  lived  (less  than 
1.5  microsecond)  evaporation  residues  are  those  of  their  respective 
daughter  functions  (4n  and  5n),  since  as  shown  in  the  ^C1  +  ^®®Tm 
reaction,  parent  and  daughter  functions  exhibit  identical  excitation 
functions  (after  correction  for  alpha  branching  ratios  less  then  1.0). 

The  tungsten  reaction  exhibits  many  anomalies  compared  to  the 
four  other  reactions  studied.  In  all  other  reactions  the  odd-even 
chlorine  beams  fused  with  odd-even  targets.  In  the  tungsten  case  the 
target  is  even-even,  however.  Utilizing  these  shell  effects,  the  ALICE 
code  has  predicted  higner  cross  sections  for  the  tungsten  reaction  in 
all  parameter  sets.  (See  Figure  54.)  Experimentally,  there  is  no 
reason  to  believe  that  the  peak  cross  section  was  obtained,  since  the 
total  evaporation  residue  cross  section  function  is  just  beginning 
Its  steep  climb.  (See  Figure  5*2.)  There  is  also  doubt  as  to  the 
identification  of  the  individual  residues.  This  has  a  drastic  effect 
on  cross  section  determination  if  a  misjudgement  is  made  between  a 
neutron  and  an  alpha  evaporation.  As  mentioned  previously,  only  a 
detailed  velocity  and/or  angular  distribution  can  differentiate  the 
two.  Due  to  time  constraints,  this  was  not  accomplished.  The 
oddities  presented,  the  1  value  dependence,  the  large  3na  evaporation, 
the  peak  cross  section  occurring  at  a  much  higher  lab  energy  as  well 
as  the  problems  of  short  half  life  for  the  parent  residues  and  the 
paucity  of  data  on  the  Protactinium  line,  make  this  a  good  candidate 


-108- 


to  pursue  further  in  heavy  ion  reaction  studies.  The  use  of  183W  as 
a  target,  although  difficult  to  obtain,  would  assist  in  determining  the 
extent  of  the  shell  effects  in  the  stability  of  the  heavy  compound 
nuclei.  While  the  data  is  insufficient  at  this  time,  the  increased 
cross  section  at  equal  Z  /A  for  odd-even  compound  nuclei  versus  even-even 
compound  systems  may  prove  significant  in  the  search  for  heavier  systems. 


Relative  Counts 
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37CI  +  181  Tq  -  2i8Th 
Angular  Distribution 


•  7.199  3np 
o  7.351  4np 
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37 Cl  +  181  To  -  2l8Th 


FIGURE  43 


FIGURE  46 
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FIGURE  48 


Relative  *  Counts 


-117- 


Angular  Distribution 


FIGURE  49 
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TABLE  V 

37  1 

Absolute  Cross  Sections  (in  yb)  for  Cl  +  Ta 


ELab(HeV> 

S0(MeV) 

155 

160 

165 

170 

175 

180 

185 

190 

4n/5n 

.23 

1.24 

1.32 

1.11 

.34 

.40 

.0 

7.655 

±, 

.04 

±.20 

±.21 

±.18 

±.05 

±.06 

±.0i 

2na/3na 

.12 

2 

.76 

1.43 

1.31 

.66 

.39 

6.899 

±.06 

±1, 

.39 

±.72 

±.66 

±.33 

±.20 

4na/5na 

.67 

1.16 

1.41 

4.98 

3.05 

4.0 

7.01 

±.22 

±.39 

±.47 

±1.66 

±1.02 

±1.3 

3np 

.42 

.80 

.78 

.77 

.07 

.0' 

7.199 

±.11 

±.20 

±.20 

±.20 

±.02 

±.0 

4np 

.14 

1.15 

2.58 

2.57 

1.8 

7.351 

±.02 

±.20 

±.45 

±.45 

±.3 

°TOT 

.12 

2 

.99 

3.76 

4.73 

5.11 

9.06 

6.09 

5.9 

±.06 

±1 

.44 

±1.25 

±1.48 

±1.38 

±2.56 

±1.55 

±1.6 
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CHAPTER  V 

RESULTS  AND  CONCLUSIONS 


2 

Figure  53  shows, the  relationship  between  the  Z  /A  value  of 

the  compound  nucleus  and  the  peak  absolute  cross  sections  of  evaporation 

residue  formation  of  the  five  reactions  (see  Table  VII).  The 

reaction  cross  sections  fall  drastically.  For  every  .5  increase  in 
2 

Z  /A  there  is  an  order  of  magnitude  drop  in  the  amount  of  compound 

nuclei  that  survive  fission.  Each  point  in  the  figure  represents  the 

peak  reaction  cross  section.  The  peak  value  is  used  since  it  was 

2 

found  that  the  relationship  between  a  and  Z  /A  is  somewhat  insensitive 
to  the  effect  of  compound  nucleus  excitation  energy  given  this  series 
of  reactions  with  similar  beams.  In  this  series  of  reactions  all  of 
the  peaks  occurred  within  a  small  bracket  of  excitation  energy.  This 
drastic  crop  in  cross  section  can  not  continue  if  more  heavy,  or 
superheavy,  nuclei  are  to  be  studied.  The  use  of  chlorine  as  a 
projectile  imparts  a  high  excitation  energy  to  the  compound  nucleus 
and  a  choice  of  an  ion  beam  at  the  Acn/4  level  may  reverse  this  trend. 
There  must  be  a  sufficient  beam  current  however,  and  these  are  not 
yet  available  in  present  systems. 

A  regression  analysis  performed  on  the  experimental 

2 

a  vs.  Z  /A  function  resulted  in  the  following  function: 

In  o  »  161.65  -  4.293  Z2/A 
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whose  "f"  value  (the  measure  of  the  reliability  of  the  hypothetical 
linear  relationship)  is  very  significant  (99.5%  reliable). 

In  the  conduct  of  the  experiments  of  this  study,  several  new 

isotopes  were  discovered.  In  the  ^C1  +  ^Tm  reaction,  the  new 
199  199n 

isotopes  Rn  and  Rn  were  found  to  have  alpha  decay  energies  of 
6.99  and  7.06  MeV,  respectively. 

In  order  to  predict  the  best  combination  of  projectile  and 
target  nucleus  to  use  in  forming  a  compound  system  that  survives 
fission  to  a  measurable  extent,  a  satisfactory  set  of  parameters  must 
be  obtained  for  statistical  evaporation  codes.  The  code  ALICE  has 
been  applied  to  the  five  reactions  in  this  study  to  determine  a  useable 
parameter  set  in  order  to  predict  future  reactions.  Several  parameter 
sets  were  tried  at  a  constant  excitation  of  the  compound  nucleus  and 
these  appear  in  figure  54  and  Table  I.  The  parameter  CLO  (the  ratio 
of  single  particle  level  densities  afission/aneutron)  1S  by  far  the 
most  sensitive  of  all  ALICE  inputs.  The  best  fit  found  was  with 
CLD  set  to  1.03.  This  continues  the  trend  announced  by  Beckerman  and 
Blann  (see  Table  VIII)  that  determined  the  adequacy  of  CLD  *  1.03  for 
lighter  systems.  This  particular  ALICE  study  concentrated  only  on 
the  total  evaporation  residue  cross  section.  Individual  evaporation 
channels  do  not  match  experimental  values  well.  Neutron  evaporations 
are  very  supressed  and  alpha  evaporations  do  not  occur  at  appropriate 
lab  energies.  Xnp  evaporation  modes  are  accurate  however.  Also,  the 
ALICE  code  predicts  several  modes  that  are  seldom  seen  in  heavy  ion 
reactions  (i.e.  3pln).  ALICE  can  also  be  used  to  generate  the  i-crit 
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TABLE  VII 

2 

Peak  Evaporation  Residue  Cross  Sections  vs.  Z  /A 


REACTION 

i!za 

CROSS  SECTION  (ub) 

37C1  +  169Tm  -  206Rn 

35.90 

1600  ±  300 

35C1  +  169Tm  -  204Rn 

36.25 

487  t  109 

37C1  +  175Lu  -*•  212Ra 

36.50 

139  ±  27 

37  Cl  +  181Ta  218Th 

37.16 

9.06  ±  2.55 

37 Cl  +  186W  ♦  223Pa 

37.13 

8.04  ±  44 

FIGURE  53 


TABLE  VIII 


*From  Beckerman  and  Blann,  PRC  May  1978,  p.  1615 
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values  for  use  in  more  sophisticated  statistical  evaporation  codes 

such  as  CASCADE  or  MSI I.  The  £-crit  values  change  for  different  lab 

energies  and  figure  5  shows  the  total  dcr/dl  versus  l  plots  for  the 

reactions  studied  only  at  the  peak  experimental  cross  section  value. 

The  corresponding  portions  of  the  total  fusion  cross  section  that 

survive  fission  are  shown  in  figure  6.  Note  the  scale  changes  and  the 

multiplicative  factor  of  the  evaporation  portions.  This  demonstrates 
-5  -4 

that  only  about  10  to  10  of  the  systems  survives  fission 
depending  on  the  reaction  studied.  It  also  shows  that  at  high 
angular  momentum  states,  here  with  an  L  value  greater  than  30,  the 
chance  of  a  compound  nucleus  surviving  fission  is  negligible.  Below 
30,  it  demonstrates  that  fission  is  still  the  dominant  process  at 
all  values  with  evaporation  remaining  only  a  very  small  part  of  the 
total  cross  section. 

ALICE,  then,  has  been  shown  to  be  an  effective  tool  in  the 
prediction  of  total  evaporation  residue  formation  in  heavy  ion  fusion 
reactions.  It  is  not  adequate  in  the  prediction  of  specific  modes  of 
evaporation,  however.  Parameters  can  be  selected  that  match  various 
specific  evaporations,  but  there  has  been  little  success  in  matching 
total  and  specific  evaporations. 

The  purpose  of  this  study  has  been  to  chart  the  expected 

2 

decrease  in  evaporation  cross  section  with  an  increase  in  Z  /A.  The 
proton  number  was  varied  from  86  to  91  and  the  atomic  number  from 
204  to  223.  In  this  small  range  from  Z^/A  of  35.9  to  37.2,  the  cross 
section  for  the  survival  against  fission  of  a  compound  nucleus 


decreased  by  a  factor  of  200.  The  ALICE  code,  using  10  different 
parameter  sets,  correctly  paralleled  this  decrease  and  confirmed  in 
the  heavier  region,  a  parameter  set  used  with  success  in  the  lighter 
systems.  The  trends  charted  in  this  study  must  be  reversed  or  the 

2 

detection  machinery  improved  if  heavier  systems  with  even  higher  Z  /A 
are  to  be  studied.  Currently  the  Heavy  Ion  Group  has  initiated  projects 
to  determine  the  dependence  on  excitation  energy  (cold  nucleus)  and 
on  shell  effects  (40Ca  +  160»162>'l63>164oy) t  since  there  seems  to  be 
a  dependence  upon  both  the  configuration  and  the  excitation  in  regards 
to  the  stability  of  the  compound  system.  Hopefully  this  study  has 
shown  the  feasibility  of  continued  study  of  the  production  of  heavy 
compound  systems  and  the  possibility  of  continuation  of  that  study 
into  heavier  domains. 
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